Evaluation of the availability of single-position treatment with a rotating gantry and the validity of deformable image registration dose assessment for pancreatic cancer carbon-ion radiotherapy.

BACKGROUND: This study aimed to evaluate the clinical acceptability of rotational gantry-based single-position carbon-ion radiotherapy (CIRT) to reduce the gastrointestinal (GI) dose in pancreatic cancer. We also evaluated the usefulness of the deformable image registration (DIR)-based dosimetry method for CIRT. MATERIAL AND METHODS: Fifteen patients with pancreatic cancer were analyzed. The treatment plans were developed for four beam angles in the supine (SP plan) and prone (PR plan) positions. In the case of using multiple positions, the treatment plan was created with two angles for each of the supine and prone position (SP + PR plan). Dose evaluation for multiple positions was performed in two ways: by directly adding the values of the DVH parameters for each position treatment plan (DVH sum), and by calculating the DVH parameters from the accumulative dose distribution created using DIR (DIR sum). The D2cc and D6cc of the stomach and duodenum were recorded for each treatment plan and dosimetry method and compared. RESULTS: There were no significant differences among any of the treatment planning and dosimetry methods (p > 0.05). The DVH parameters for the stomach and duodenum were higher in the PR plan and SP plan, respectively, and DVH sum tended to be between the SP and PR plans. DVH sum and DIR sum, DVH sum tended to be higher for D2cc and DIR sum tended to be higher for D6cc . CONCLUSION: There were no significant differences in the GI dose, which suggests that treatment with a simple workflow performed in one position should be clinically acceptable. In CIRT, DIR-based dosimetry should be carefully considered because of the potential for increased uncertainty due to the steep dose distributions.

A robust treatment planning approach for chest motion in postmastectomy chest wall intensity modulated radiation therapy.

PURPOSE: Chest wall postmastectomy radiation therapy (PMRT) should consider the effects of chest wall respiratory motion. The purpose of this study is to evaluate the effectiveness of robustness planning intensity modulated radiation therapy (IMRT) for respiratory movement, considering respiratory motion as a setup error. MATERIAL AND METHODS: This study analyzed 20 patients who underwent PMRT (10 left and 10 right chest walls). The following three treatment plans were created for each case and compared. The treatment plans are a planning target volume (PTV) plan (PP) that covers the PTV within the body contour with the prescribed dose, a virtual bolus plan (VP) that sets a virtual bolus in contact with the body surface and prescribing the dose that includes the PTV outside the body contour, and a robust plan (RP) that considers respiratory movement as a setup uncertainty and performs robust optimization. The isocenter was shifted to reproduce the chest wall motion pattern and the doses were recalculated for comparison for each treatment plan. RESULT: No significant difference was found between the PP and the RP in terms of the tumor dose in the treatment plan. In contrast, VP had 3.5% higher PTV Dmax and 5.5% lower PTV V95% than RP (p < 0.001). The RP demonstrated significantly higher lung V20Gy and Dmean by 1.4% and 0.4 Gy, respectively, than the PP. The RP showed smaller changes in dose distribution affected by chest wall motion and significantly higher tumor dose coverage than the PP and VP. CONCLUSION: We revealed that the RP demonstrated comparable tumor doses to the PP in treatment planning and was robust for respiratory motion compared to both the PP and the VP. However, the organ at risk dose in the RP was slightly higher; therefore, its clinical use should be carefully considered.

Error on the stopping power ratio of ERKODENT's mouthpiece for head and neck carbon ion radiotherapy treatment.

The errors on the stopping power ratio (SPR) of mouthpiece samples from ERKODENT were evaluated. Erkoflex and Erkoloc-pro from ERKODENT and samples that combined Erkoflex and Erkoloc-pro were computed tomography (CT)-scanned using head and neck (HN) protocol at the East Japan Heavy Ion Center (EJHIC), and the values were averaged to obtain the CT number. The integral depth dose of the Bragg curve with and without these samples was measured for 292.1, 180.9, and 118.8 MeV/u of the carbon-ion pencil beam using an ionization chamber with concentric electrodes at the horizontal port of the EJHIC. The average value of the water equivalent length (WEL) of each sample was obtained from the difference between the range of the Bragg curve and the thickness of the sample. To calculate the difference between the theoretical and measured values, the theoretical CT number and SPR value of the sample were calculated using the stoichiometric calibration method. Compared with the Hounsfield unit (HU)-SPR calibration curve used at the EJHIC, the SPR error on each measured and theoretical value was calculated. The WEL value of the mouthpiece sample had an error of approximately 3.5% in the HU-SPR calibration curve. From this error, it was evaluated that for a mouthpiece with a thickness of 10 mm, a beam range error of approximately 0.4 mm can occur, and for a mouthpiece with a thickness of 30 mm, a beam range error of approximately 1 mm can occur. For a beam passing through the mouthpiece in HN treatment, it would be practical to consider a mouthpiece margin of 1 mm to avoid beam range errors if ions pass through the mouthpiece.

Robust treatment planning in scanned carbon-ion radiotherapy for pancreatic cancer: Clinical verification using in-room computed tomography images.

Purpose: Carbon-ion beam (C-beam) has a sharp dose distribution called the Bragg peak. Carbon-ion radiation therapy, such as stereotactic body radiotherapy in photon radiotherapy, can be completed in a short period by concentrating the radiation dose on the tumor while minimizing the dose to organs at-risk. However, the stopping position of C-beam is sensitive to density variations along the beam path and such variations can lower the tumor dose as well as cause the delivery of an unexpectedly high dose to the organs at risk. We evaluated the clinical efficacy of a robust planning technique considering gastrointestinal gas (G-gas) to deliver accurate radiation doses in carbon-ion radiotherapy for pancreatic cancer. Materials and methods: We focused on the computed tomography (CT) value replacement method. Replacement signifies the overwriting of CT values in the CT images. The most effective replacement method for robust treatment planning was determined by verifying the effects of the three replacement patterns. We selected 10 consecutive patients. Pattern 1 replaces the CT value of the G-gas contours with the value of the region without G-gas (P1). This condition indicates a no-gas state. Pattern 2 replaces each gastrointestinal contour using the mean CT value of each contour (P2). The effect of G-gas was included in the replacement value. Pattern 3 indicates no replacement (P3). We analyzed variations in the target coverage (TC) and homogeneity index (HI) from the initial plan using in-room CT images. We then performed correlation analysis on the variations in G-gas, TC, and HI to evaluate the robustness against G-gas. Results: Analysis of variations in TC and HI revealed a significant difference between P1 and P3 and between P2 and P3. Although no statistically significant difference was observed between P1 and P2, variations, including the median, tended to be fewer in P2. The correlation analyses for G-gas, TC, and HI showed that P2 was less likely to be affected by G-gas. Conclusion: For a treatment plan that is robust to G-gas, P2 mean replacement method should be used. This method does not necessitate any particular software or equipment, and is convenient to implement in clinical practice.

Development of a quantitative analysis method for assessing patient body surface deformation using an optical surface tracking system.

This study aimed to develop a new method to quantitatively analyze body shape changes in patients during radiotherapy without additional radiation exposure using an optical surface tracking system. This method's accuracy was evaluated using a cubic phantom with a known shift. Surface images of three-dimensionally printed phantoms, which simulated the head and neck shapes of real patients before and after treatment, were used to create a deformation surface area histogram. The near-maximum deformation value covering an area of 2 cm2 in the surface image (Def-2cm2) was calculated. A volumetric modulated arc therapy (VMAT) plan was also created on the pre-treatment phantom, and the dose distribution was recalculated on the post-treatment phantom to compare the dose indices. Surface images of four patients were analyzed to evaluate Def-2cm2 and examine whether this method can be used in clinical cases. Experiments with the cubic phantom resulted in a mean deformation error of 0.08 mm. With head and neck phantoms, the Def-2cm2 value was 17.5 mm, and the dose that covered 95% of the planning target volume in the VMAT plan decreased by 11.7%, indicating that deformation of the body surface may affect the dose distribution. Although analysis of the clinical data showed no clinically relevant deformation in any of the cases, slight skin sagging and respiratory changes in the body surface were observed. The proposed method can quantitatively and accurately evaluate the deformation of a body surface. This method is expected to be used to make decisions regarding modifications to treatment plans.

Prediction of the minimum spacer thickness required for definitive radiotherapy with carbon ions and photons for pelvic tumors: an in silico planning study using virtual spacers.

We aimed to predict the minimum distance between a tumor and the gastrointestinal (GI) tract that can satisfy the dose constraint by creating simulation plans with carbon-ion (C-ion) radiotherapy (RT) and photon RT for each case assuming insertion of virtual spacers of various thicknesses. We enrolled 55 patients with a pelvic tumor adjacent to the GI tract. Virtual spacers were defined as the overlap volume between the GI tract and the volume expanded 7-17 mm from the gross tumor volume (GTV). Simulation plans (70 Gy in 35 fractions for at least 95% of the planning target volume [PTV]) were created with the lowest possible dose to the GI tract under conditions that meet the dose constraints of the PTV. We defined the minimum thickness of virtual spacers meeting D2 cc of the GI tract <50 Gy as 'MTS'. Multiple regression was used with explanatory variables to develop a model to predict MTS. We discovered that MTSs were at most 9 mm and 13 mm for C-ion RT and photon RT plans, respectively. The volume of overlap between the GI tract and a virtual spacer of 14 mm in thickness (OV14)-PTV was found to be the most important explanatory variable in the MTS prediction equation for both C-ion and photon RT plans. Multiple R2 values for the regression model were 0.571 and 0.347 for C-ion RT and photon RT plans, respectively. In conclusion, regression equations were developed to predict MTS in C-ion RT and photon RT.

Comparison of predictive performance for toxicity by accumulative dose of DVH parameter addition and DIR addition for cervical cancer patients.

We compared predictive performance between dose volume histogram (DVH) parameter addition and deformable image registration (DIR) addition for gastrointestinal (GI) toxicity in cervical cancer patients. A total of 59 patients receiving brachytherapy and external beam radiotherapy were analyzed retrospectively. The accumulative dose was calculated by three methods: conventional DVH parameter addition, full DIR addition and partial DIR addition. ${D}_{2{cm}^3}$, ${D}_{1{cm}^3}$ and ${D}_{0.1{cm}^3}$ (minimum doses to the most exposed 2 cm3, 1cm3 and 0.1 cm3 of tissue, respectively) of the rectum and sigmoid were calculated by each method. V50, V60 and V70 Gy (volume irradiated over 50, 60 and 70 Gy, respectively) were calculated in full DIR addition. The DVH parameters were compared between toxicity (≥grade1) and non-toxicity groups. The area under the curve (AUC) of the receiver operating characteristic (ROC) curves were compared to evaluate the predictive performance of each method. The differences between toxicity and non-toxicity groups in ${D}_{2{cm}^3}$ were 0.2, 5.7 and 3.1 Gy for the DVH parameter addition, full DIR addition and partial DIR addition, respectively. The AUCs of ${D}_{2{cm}^3}$ were 0.51, 0.67 and 0.57 for DVH parameter addition, full DIR addition and partial DIR addition, respectively. In full DIR addition, the difference in dose between toxicity and non-toxicity was the largest and AUC was the highest. AUCs of V50, V60 and V70 Gy were 0.51, 0.63 and 0.62, respectively, and V60 and V70 were high values close to the value of ${D}_{2{cm}^3}$ of the full DIR addition. Our results suggested that the full DIR addition may have the potential to predict toxicity more accurately than the conventional DVH parameter addition, and that it could be more effective to accumulate to all pelvic irradiation by DIR.

The impact of 4DCT-ventilation imaging-guided proton therapy on stereotactic body radiotherapy for lung cancer.

Functional lung avoidance during radiotherapy can help reduce pulmonary toxicity. This study assessed the potential impact of four-dimensional computed tomography (4DCT)-ventilation imaging-guided proton radiotherapy (PT) on stereotactic body radiotherapy (SBRT) by comparing it with three-dimensional conformal radiotherapy (3D-CRT) and volumetric modulated arc therapy (VMAT), which employ photon beams. Thirteen lung cancer patients who received SBRT with 3D-CRT were included in the study. 4DCT ventilation was calculated using the patients' 4DCT data, deformable image registration, and a density-change-based algorithm. Three functional treatment plans sparing the functional lung regions were developed for each patient using 3D-CRT, VMAT, and PT. The prescribed doses and dose constraints were based on the Radiation Therapy Oncology Group 0618 protocol. We evaluated the region of interest (ROI) and functional map-based dose-function metrics for 4DCT ventilation and the irradiated dose. Using 3D-CRT, VMAT, and PT, the percentages of the functional lung regions that received ≥ 5 Gy (fV5) were 26.0%, 21.9%, and 10.7%, respectively; the fV10 were 14.4%, 11.4%, and 9.0%, respectively; and fV20 were 6.5%, 6.4%, and 6.6%, respectively, and the functional mean lung doses (fMLD) were 5.6 Gy, 5.2 Gy, and 3.8 Gy, respectively. These results indicated that PT resulted in a significant reduction in fMLD, fV5, and fV10, but not fV20. The use of PT reduced the radiation to highly functional lung regions compared with those for 3D-CRT and VMAT while meeting all dose constraints.

Quantitative analysis of intra-fractional variation in CT-based image guided brachytherapy for cervical cancer patients.

We quantified intra-fractional dose variation and organ movement during CT-based 3D-image guided brachytherapy (3D-IGBT) in cervical cancer patients. Fifteen patients who underwent CT-based 3D-IGBT were studied. For all patients, pre-delivery CT for treatment planning after applicator insertion and post-delivery CT after dose delivery without changing the applicator position were acquired. Pre- and post-delivery CT were rigidly fused by matching the inserted applicator and planned dose on pre-delivery CT (pre-delivery dose) was mapped on post-delivery CT (post-delivery dose). D2, D1, and D0.1 cm3 of the rectum and bladder were compared between pre- and post-delivery doses with contours on each CT image. Organ movement and deformation was evaluated using deformation vector fields calculated by deformable image registration between pre- and post-delivery CT. We also evaluated dose variation and DVF between with and without a catheter to control filling. Differences in all DVH parameters were <±3% in physical dose and ± 5% in EQD2. However, a > 15% dose difference was found in 13.8% of the fractions in rectum D2 cm3 and in 11.1% of those in bladder D2 cm3. The mean value of DVF for bladder was larger than that of rectum, especially for the superior-inferior (S-I) direction. Insertion catheters in bladder reduced mean dose and DVF variation compared with that of without catheters. In fraction groups with large dose increasing, DVF in the S-I direction was significantly larger than that of other fraction groups. Our results indicated that preparation is needed to reduce changes in the S-I direction affect dose variation.

In silico comparison of the dosimetric impacts of a greater omentum spacer for abdominal and pelvic tumors in carbon-ion, proton and photon radiotherapy.

PURPOSE: The purpose of this study was to compare carbon-ion (C-ion), proton and photon radiotherapy (RT) plans with regard to dose reduction of the gastrointestinal (GI) tract by using a greater omentum spacer (GO spacer). METHODS: We retrospectively retrieved data for ten patients who received the GO spacer as surgical spacer placement for abdominal and pelvic tumors. Simulation plans were created on pre-spacer Computed Tomography (CT) and post-spacer CT for C-ion RT, proton RT and photon RT to compare the dose of the GI tract. The plans were normalized so that at least 95% of the planning target volume (PTV) received 70 Gy (relative biological effectiveness equivalent) delivered in 35 fractions. All plans were created with the lowest possible dose to the GI tract under conditions that meet the dose constraints for the PTV and spinal cord (maximum dose < 45 Gy). The part of the GI tract to be evaluated was defined as that most adjacent to the PTV. C-ion RT plans and proton RT plans were calculated by a spot scanning technique, and photon RT plans were calculated employing by fixed-field intensity-modulated radiation therapy. RESULTS: D2 cc and V10-70 of the GI tract were significantly lower on post-spacer plans than on pre-spacer plans for all three RT modalities. Regarding post-spacer plans, D2 cc of the GI tract was significantly lower on C-ion RT plans and proton RT plans than on photon RT plans (C-ion vs photon p = 0.001, proton vs photon p = 0.002). However, there was no significant difference between C-ion RT plans and proton RT plans for D2 cc of the GI tract (C-ion vs proton p = 0.992). In the photon RT plan for one patient, D2 cc of the GI tract did not meet < 50 Gy. CONCLUSIONS: The GO spacer shows a significant dose reduction effect on the GI tract.

Evaluation of functionally weighted dose-volume parameters for thoracic stereotactic ablative radiotherapy (SABR) using CT ventilation.

For the purpose of reducing radiation pneumontisis (RP), four-dimensional CT (4DCT)-based ventilation can be used to reduce functionally weighted lung dose. This study aimed to evaluate the functionally weighted dose-volume parameters and to investigate an optimal weighting method to realize effective planning optimization in thoracic stereotactic ablative radiotherapy (SABR). Forty patients treated with SABR were analyzed. Ventilation images were obtained from 4DCT using deformable registration and Hounsfield unit-based calculation. Functionally-weighted mean lung dose (fMLD) and functional lung fraction receiving at least x Gy (fVx) were calculated by two weighting methods: thresholding and linear weighting. Various ventilation thresholds (5th-95th, every 5th percentile) were tested. The predictive accuracy for CTCAE grade ≧ 2 pneumonitis was evaluated by area under the curve (AUC) of receiver operating characteristic analysis. AUC values varied from 0.459 to 0.570 in accordance with threshold and dose-volume metrics. A combination of 25th percentile threshold and fV30 showed the best result (AUC: 0.570). AUC values with fMLD, fV10, fV20, and fV40 were 0.541, 0.487, 0.548 and 0.563 using a 25th percentile threshold. Although conventional MLD, V10, V20, V30 and V40 showed lower AUC values (0.516, 0.477, 0.534, 0.552 and 0.527), the differences were not statistically significant. fV30 with 25th percentile threshold was the best predictor of RP. Our results suggested that the appropriate weighting should be used for better treatment outcomes in thoracic SABR.

Impact of a commercially available model-based dose calculation algorithm on treatment planning of high-dose-rate brachytherapy in patients with cervical cancer.

We evaluated the impact of model-based dose calculation algorithms (MBDCAs) on high-dose-rate brachytherapy (HDR-BT) treatment planning for patients with cervical cancer. Seven patients with cervical cancer treated using HDR-BT were studied. Tandem and ovoid applicators were used in four patients, a vaginal cylinder in one, and interstitial needles in the remaining two patients. MBDCAs were applied to the Advanced Collapsed cone Engine (ACE; Elekta, Stockholm, Sweden). All plans, which were originally calculated using TG-43, were re-calculated using both ACE and Monte Carlo (MC) simulations. Air was used as the rectal material. The mean difference in the rectum D2cm3 between ACErec-air and MCrec-air was 8.60 ± 4.64%, whereas that in the bladder D2cm3 was -2.80 ± 1.21%. Conversely, in the small group analysis (n = 4) using water instead of air as the rectal material, the mean difference in the rectum D2cm3 between TG-43 and ACErec-air was 11.87 ± 2.65%, whereas that between TG-43 and ACErec-water was 0.81 ± 2.04%, indicating that the use of water as the rectal material reduced the difference in D2cm3 between TG-43 and ACE. Our results suggested that the differences in the dose-volume histogram (DVH) parameters of TG-43 and ACE were large for the rectum when considerable air (gas) volume was present in it, and that this difference was reduced when the air (gas) volume was reduced. Also, ACE exhibited better dose calculation accuracy than that of TG-43 in this situation. Thus, ACE may be able to calculate the dose more accurately than TG-43 for HDR-BT in treating cervical cancers, particularly for patients with considerable air (gas) volume in the rectum.

Evaluating the Toxicity Reduction With Computed Tomographic Ventilation Functional Avoidance Radiation Therapy.

PURPOSE: Computed tomographic (CT) ventilation imaging is a new modality that uses 4-dimensional (4D) CT information to calculate lung ventilation. Although retrospective studies have reported on the reduction in dose to functional lung, no work to our knowledge has been published in which the dosimetric improvements have been translated to a reduction in the probability of pulmonary toxicity. Our work estimates the reduction in toxicity for CT ventilation-based functional avoidance planning. METHODS AND MATERIALS: Seventy previously treated lung cancer patients who underwent 4DCT imaging were used for the study. CT ventilation maps were calculated with 4DCT deformable image registration and a density change-based algorithm. Pneumonitis was graded on the basis of imaging and clinical presentation. Maximum likelihood methods were used to generate normal tissue complication probability (NTCP) models predicting grade 2 or higher (2+) and grade 3+ pneumonitis as a function of dose (V5 Gy, V10 Gy, V20 Gy, V30 Gy, and mean dose) to functional lung. For 30 patients a functional plan was generated with the goal of reducing dose to the functional lung while meeting Radiation Therapy Oncology Group 0617 constraints. The NTCP models were applied to the functional plans and the clinically used plans to calculate toxicity reduction. RESULTS: By the use of functional avoidance planning, absolute reductions in grade 2+ NTCP of 6.3%, 7.8%, and 4.8% were achieved based on the mean fV20 Gy, fV30 Gy, and mean dose to functional lung metrics, respectively. Absolute grade 3+ NTCP reductions of 3.6%, 4.8%, and 2.4% were achieved with fV20 Gy, fV30 Gy, and mean dose to functional lung. Maximum absolute reductions of 52.3% and 16.4% were seen for grade 2+ and grade 3+ pneumonitis for individual patients. CONCLUSION: Our study quantifies the possible toxicity reduction from CT ventilation-based functional avoidance planning. Reductions in grades 2+ and 3+ pneumonitis were 7.1% and 4.7% based on mean dose-function metrics, with reductions as high as 52.3% for individual patients. Our work provides seminal data for determining the potential toxicity benefit from incorporating CT ventilation into thoracic treatment planning.

Evaluation of rectum and bladder dose accumulation from external beam radiotherapy and brachytherapy for cervical cancer using two different deformable image registration techniques.

We evaluated dose-volume histogram (DVH) parameters based on deformable image registration (DIR) between brachytherapy (BT) and external beam radiotherapy (EBRT) that included a center-shielded (CS) plan. Eleven cervical cancer patients were treated with BT, and their pelvic and CS EBRT were studied. Planning CT images for EBRT and BT (except for the first BT, used as the reference image) were deformed with DIR to reference image. We used two DIR parameter settings: intensity-based and hybrid. Mean Dice similarity coefficients (DSCs) comparing EBRT with the reference for the uterus, rectum and bladder were 0.81, 0.77 and 0.83, respectively, for hybrid DIR and 0.47, 0.37 and 0.42, respectively, for intensity-based DIR (P < 0.05). D1 cm3 for hybrid DIR, intensity-based DIR and DVH addition were 75.1, 81.2 and 78.2 Gy, respectively, for the rectum, whereas they were 93.5, 92.3 and 94.3 Gy, respectively, for the bladder. D2 cm3 for hybrid DIR, intensity-based DIR and DVH addition were 70.1, 74.0 and 71.4 Gy, respectively, for the rectum, whereas they were 85.4, 82.8 and 85.4 Gy, respectively, for the bladder. Overall, hybrid DIR obtained higher DSCs than intensity-based DIR, and there were moderate differences in DVH parameters between the two DIR methods, although the results varied among patients. DIR is only experimental, and extra care should be taken when comparing DIR-based dose values with dose-effect curves established using DVH addition. Also, a true evaluation of DIR-based dose accumulation would require ground truth data (e.g. measurement with physical phantom).

Evaluation of deformable image registration between external beam radiotherapy and HDR brachytherapy for cervical cancer with a 3D-printed deformable pelvis phantom.

PURPOSE: In this study, we developed a 3D-printed deformable pelvis phantom for evaluating spatial DIR accuracy. We then evaluated the spatial DIR accuracies of various DIR settings for cervical cancer. METHODS: A deformable female pelvis phantom was created based on patient CT data using 3D printing. To create the deformable uterus phantom, we first 3D printed both a model of uterus and a model of the internal cavities of the vagina and uterus. We then made a mold using the 3D printed uterus phantom. Finally, urethane was poured into the mold with the model of the internal cavities in place, creating the deformable uterus phantom with a cavity into which an applicator could be inserted. To create the deformable bladder phantom, we first 3D printed models of the bladder and of the same bladder scaled down by 2 mm. We then made a mold using the larger bladder model. Finally, silicone was poured into the mold with the smaller bladder model in place to create the deformable bladder phantom with a wall thickness of 2 mm. To emulate the anatomical bladder, water was poured into the created bladder. We acquired phantom image without applicator for EBRT. Then, we inserted the applicator into the phantom to simulate BT. In this situation, we scanned the phantom again to obtain the phantom image for BT. We performed DIR using the two phantom images in two cases: Case A, with full bladder (170 ml) in both EBRT and BT images; and Case B with full bladder in the BT image and half-full bladder (100 ml) in the EBRT image. DIR was evaluated using Dice similarity coefficients (DSCs) and 31 landmarks for the uterus and 25 landmarks for the bladder. A hybrid intensity and structure DIR algorithm implemented in RayStation with four DIR settings was evaluated. RESULTS: On visual inspection, reasonable agreement in shape of the uterus between the phantom and patient CT images was observed for both EBRT and BT, although some regional disagreements in shape of the bladder and rectum were apparent. The created phantom could reproduce the actual patient's uterus deformation by the applicator. For both Case A and B, large variation was seen in landmark error among the four DIR parameters. In addition, although DSCs were comparable, moderate differences in landmark error existed between the two different DIR parameters selected from the four DIR parameters (i.e., DSC = 0.96, landmark error = 13.2 ± 5.7 mm vs. DSC = 0.97, landmark error = 9.7 ± 4.0 mm). This result suggests that landmark error evaluation might thus be more effective than DSC for evaluating DIR accuracy. CONCLUSIONS: Our developed phantom enabled the evaluation of spatial DIR accuracy for the female pelvic region for the first time. Although the DSCs are high, the spatial errors can still be significant and our developed phantom facilitates their quantification. Our results showed that optimization is needed to identify suitable DIR settings. For determining suitable DIR settings, our method of evaluating spatial DIR accuracy using the 3D-printed phantom may prove helpful.

Evaluation of the performance of deformable image registration between planning CT and CBCT images for the pelvic region: comparison between hybrid and intensity-based DIR.

This study aimed to evaluate the performance of the hybrid deformable image registration (DIR) method in comparison with intensity-based DIR for pelvic cone-beam computed tomography (CBCT) images, using intensity and anatomical information. Ten prostate cancer patients treated with intensity-modulated radiation therapy (IMRT) were studied. Nine or ten CBCT scans were performed for each patient. First, rigid registration was performed between the planning CT and all CBCT images using gold fiducial markers, and then DIR was performed. The Dice similarity coefficient (DSC) and center of mass (COM) displacement were used to evaluate the quantitative DIR accuracy. The average DSCs for intensity-based DIR for the prostate, rectum, bladder, and seminal vesicles were 0.84 ± 0.05, 0.75 ± 0.05, 0.69 ± 0.07 and 0.65 ± 0.11, respectively, whereas those values for hybrid DIR were 0.98 ± 0.00, 0.97 ± 0.01, 0.98 ± 0.00 and 0.94 ± 0.03, respectively (P < 0.05). The average COM displacements for intensity-based DIR for the prostate, rectum, bladder, and seminal vesicles were 2.0 ± 1.5, 3.7 ± 1.4, 7.8 ± 2.2 and 3.6 ± 1.2 mm, whereas those values for hybrid DIR were 0.1 ± 0.0, 0.3 ± 0.2, 0.2 ± 0.1 and 0.6 ± 0.6 mm, respectively (P < 0.05). These results showed that the DSC for hybrid DIR had a higher DSC value and smaller COM displacement for all structures and all patients, compared with intensity-based DIR. Thus, the accumulative dose based on hybrid DIR might be trusted as a high-precision dose estimation method that takes into account organ movement during treatment radiotherapy.